Generating the mass gap of the sine-Gordon model 
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We discuss in this study the possibility of finding a finite mass gap in the broken phase of the 
sine-Gordon model in d = 2 using the functional flows. We demonstrate that the signal of the 
presence of massive excitations, a finite positively-curved blocked potential around its minima, is 
recovered only in our treatment. The usual results based on the flow of the Fourier expansion of 
the blocked action are then shown to actually fit a singularity. 
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I. Introduction - The sine-Gordon Model (sGM) 
whose bare Lagrangian in d euclidean dimensions reads : 



(1) 



has attracted a considerable interest these last decades 
due to its unique features. In d = 2, the Coleman point 
(u c , z c ) = (0, gjjp) plays a special role in the weak-coupling 
ub limit and provides a three-sectors phase diagram : a 
massless sector for zb < z c due to a vanishing inter- 
action, and a strong-coupling one for zb > z c . These 
two sectors are separated by a cross-over regime around 
zb — z c for finite ub and the model is known to undergo 
a Kosterlitz-Thouless phase transition [1]. The sponta- 
neously broken periodicity symmetry phase corresponds 
to the two last sectors and gives rise to dynamical mass 
generation. It exhibits solitons and bound-states and is 
quantum S-dual to the massive Thirring model where its 
solitons correspond to the fermions of the latter [2] . The 
exact spectrum of these topological excitations is known 
[3]. At the partition function level, it is also dual to the 
classical neutral Coulomb gas in all dimensions and is 
thus of great interest to study the roughening transition 
[4] and the universal jump of the superfluid density in 
d = 2. Finally, the sGM is also used to describe the 
thermal fluctuations of membranes in a periodic pinning 
potential (e. g. [5]). 

Because the sGM has been exactly solved, it is an 
ideal testing ground to understand better theoretical 
approaches, such as functional Renormalization Group 
(fRG). Understanding how fRG works for the sGM will 
provide us new insights to study models where the exact 
methods are beyond our abilities but fRG is immediatly 
applicable, such as a more general periodic interaction 
an/or in another dimension. 

The present Letter focuses on the signature of the dynam- 
ical mass generation, which is the central phenomenom 
occuring in the IR of the broken phase. We will demon- 
strate that the appearance of this new mass scale can be 
seen only in our approach. 

II. Functional flows - In the following, we use the 
Effective Average Action (EAA) method. Its strategy 



is to build a continuous set of functionals interpolat- 
ing smoothly between the bare action of the theory one 
whishes to solve and its effective action. These function- 
als T k -the effective average action- are demonstrated to 
follow the flow equation : 



kd k T k = -Tr 



kd k R k (p) 



(2) 



where the regulator R k (p 2 ) is a well-chosen function of 
k and p 2 [6]. Obviously, the functional equation (2) is 
in general too difficult to be solved directly, and calls for 
further assumptions. The most common approximation 
is to project the EAA on the functional ansatz : 



d 2 x ^M d ^ d ^ + Vk[(j)] (3) 



which gives now a set of two coupled flow equations for 
Vh(<p) and Z k ((f)) (see e. g. [7]). We solved this system 
with no further approximation for the first time using 
the algorithm [8] for the power-law regulator R k (p 2 ) — 

P 2 (w^j with b > 1 and we discuss in the following how 
to compare these results with the usual perturbative ap- 
proach and how to extract new physical informations. 
a. The perturbative flow. Instead of solving directly 
the flow of V k {<j)) and Z k ((f>), the perturbative approach 
projects their evolution onto their first harmonic in the 
small u k limit [9] : 
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where the tilde stands for dimensionless quantities. The 
weak-coupling expansion is thus actually not driven by u 
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but rather by zl -i/ h ■ The coefficients can be systemati- 
cally computed analytically, and two universal -regulator 
independent- features appear. The Coleman point is al- 
ways reproduced as demonstrated in [10] i.e. A\ = j- 
and the coefficient B2 is negative, providing a positive 
anomalous dimension. The set of differential equations 
(4)- (5) yields an RG invariant r that can be cxplicitely 
computed at order u\ (and u\) e. g. for b = 2 : 

and thus provides analytically the trajectories in the 
(uk, £fc)-plane of the sytem. These trajectories can be 
casted in three categories, corresponding to the three sec- 
tors : the massless phase for r < and Zk < z c , the 
strong-coupling one for r < and Zk > z c and the cross- 
over regime for r > 0. The latter is limited by the separa- 
trix t — 0, plotted on fig. 1. On the same figure, we draw 
a perturbative trajectory of the strong-coupling sector 
and for comparison a serie of trajectories coming from our 
full solution of Vk((j)) and Zk(<fi) defining the couplings 
"fc = n f-?<f>Vk(<fi)cos(<p) and z k = ^ J^%4> Z k (<j>). 
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FIG. 1: Trajectories from the full solution of (3) projected 
onto the (ilk, Zfc) plane in the strong-coupling and the cross- 
over sectors (lines), the separatrix defined by r = (+) and 
a perturbative trajectory (x) for b — 2. 

The UV flow of the full solution and the perturbative 
result are in good agreement and one can thus use the 
RG invariant (6) to extract the correlation length of the 
dual XY model, and the vortex-vortex two-points func- 
tion critical exponent in the vicinity of the Coleman point 
[9] . On the contrary, the IR flow is not driving the system 
in the close neighborhood of the separatrix, as perturba- 
tively expected. As a result, the use of the RG invariant 
r to give a rough estimate of the soliton mass [11] is in- 
correct. 

b. Breakdown of the perturbative approach. The discrep- 
ancy between the perturbative results and the full solu- 
tion has many origins. We made a small i"i /b expan- 
sion, which turns out to be large in the IR even if both 



couplings Uk and Zk are growing as k is decreased, due 
to the /3-function of Uk winning the competition with the 
one of Zk- As a result, truncating the equation (4)-(5) at 
the second (or any finite) order turns out to be illegal in 
the IR and one has to solve the first lines of (4) and (5) 
without expansion [9]. In addition, taking into account 
only the first harmonic of period 2-7T of the potential Vk 
uses the relevance classification in the UV, neglecting the 
appearance of higher period harmonics. This is certainly 
legal at the beginning of the flow, but it is less and less 
correct as the non-linearities generate these harmonics 
and as Zk grows : below the scale k for which Zk = n 2 z c , 
the terms in oc cos(n<p) in Vk will start being relevant. 
Beyond this technical issue is actually hidden the fact 
that the IR fluctuations are expected to be dominated 
by the tunneling between different equivalent minima - 
this is the symmetry broken phase- i. e. the contribution 
of the topological excitations. Dropping terms of period 
2n7r neglects a part of the contributions that the topolog- 
ical configurations of winding number |n| should see. As 
a result, using a finite order Fourier expansion of Vk{4>) is 
dangerous. We now demonstrate that it is also the case 
for Zk (<j>) if one wants to see the effect of the mass of the 
excitations. 

III. Signature of the excitations - As it has been 
emphasized in [10] already at Local Potential Approxima- 
tion (LPA), the behaviour of the dimensionful quantities 
such as the coupling Uk are not informative at all in the 
case of the sGM, because the effective potential -that 
corresponds to Vk=o{4>)- is a constant in both phases, 
being the only possibility for a periodic and convex func- 
tion of its argument. This tells us that what counts is 
how the effective potential is reached i. e. the difference 
is made on the dimensionless potential Vk- 
The minima of the bare potential (1) are located in 
v = -K + 2mr and we expect that the quantum vacua, solu- 
tions of the equation of motion for the effective action are 
sitting at the same location. The vacua are in our case 
actually built progressively from the minima of Tk, also 
located in v. In the broken phase, for k large enough, 
we feel the contributions of the topological excitations 
which tend to increase slowly V^{v) with k (minima se- 
lection) but their mass are roughly priceless in units of k. 
For scales k lower than say |, we start feeling the mass 
of some of the excitations and any attempt of reaching 
an excited state has thus a finite cost in k, accelerat- 
ing the flow of VI' (v) (massive cross-over). Finally, for 
k « one fully feels the mass of the excitations and 
the energy needed to feed excited states is constant i. c. 
V"(v) diverges as fc -2 as k is descreased. 
a. The Fourier expansion ofVk and Zk- As a result, one 
realizes that the signature of the excitations can be seen 
only if one can carefully study the close neighborhood of 
v of Tk- Because the bare theory (1) is periodic and this 
periodicity is preserved by (2) and (3), it would seem nat- 
ural to represent Vk and Zk with their Fourier series. The 
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problem is that already at LPA, this strategy is known 
to fail for V k ((f) [10]. The reason is that the expansion 
does not converge due to the severe competition between 
the convex regions -i.e. around the minima, where the 
potential is well-behaved- and the concave ones -around 
the extrema, where the potential lives near a singularity 
see below. This abrupt change of behavior between the 
convex and the concave regions for the potential V k {4>) is 
thus expected to propagate to Z k ((f>) when we solve (3) 
that couples the two evolutions. 
b. Instability driven Fourier flow. Because the Fourier 
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FIG. 2: Fhe same trajectories as in fig. 2 projected onto the 
(— Vfc'(0), Zfc(0)) plane for b — 2. The extra trajectory traces 
the limit of validity of the approach given by a = 0, see (7). 

expansion is probably breaking down, it is informative to 
compare the value of V k (</>) and Z k (<j>) from our full solu- 
tion at various <fi and the value of the Fourier coefficients 
u k and z k one can build from them. At the special point 
<p = that places us in the middle of the concave regions, 
one can plot the trajectories in the (—V"(0), Z k (0)) plane 
[15] on fig 2 and realize that it has precisely the same 
shape as the usual phase diagram of fig. 1. The IR flow 
now gets close to the line of equation 

a = 1 + V k "(Q)Z k (Q)- 1+1 / b (b - if-^b- 1 > (7) 

(2) 

which is the limit of validity V. + Rk = for (2) under 
the assumption (3) when using our regulator. Therefore, 
the Fourier expansion of T k , even when valid, is actually 
not only mimicking the concave region, but also fitting a 
singularity and thus is not physically relevant. A simple 
explanation is provided by the anomalous dimension 

%(0) = -kd k Log(Z k (cb)) (8) 

computed for (f> = and <fi — v, and compare them with 
the one from the reconstructed Fourier coefficient z k . The 
figure 3 shows that in the UV, the behavior is the same in 
convex and concav region and in good agreement with the 
Fourier representation, but dramatically differ in the IR. 
In particular, it shows that the anomalous dimension goes 
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FIG. 3: The running of the anomalous dimensions r)k{4>) in 
= and cj> = v (lines) together with the Fourier recon- 
structed values (crosses) as a function of k for b = 2 . 



to around the minima while it tends to a large value 
in the concave regions. This large value in the concave 
regions makes them count more and more when going 
IR, which explains why they are the ones seen by the 
Fourier expansion. It also explains why even if one solves 
fully V k (</>) together with z k defined as in the first line of 
(5), one can not stabilize the flow around the minima : 
because z k - that fits the concave regions behavior- keeps 
running in the flow equation for V kl as in [5, 9]. 
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FIG. 4: The physical phase diagram : the same trajectories 
as in fig. 1 and 2 projected in the (V k (v), Z k (v)) plane for 
b — 2. The focalisation effect in (SuZkiv))- 1 ~ 0.35 is not 
genuine, as one can pick different bare theories converging to 
other IR values. 

IV. The physics around the minima - We can fi- 
nally turn our attention to see how behaves near v 
to get the physical phase diagram that is plotted in the 
(V k (v),Z k (v))-plaiie [14]. The figure 4 exhibits actually 
a similar behaviour in a much more flexible approach as 
the one obtained in a different context, using Wegner's 
blocking [12]. As expected, the values of Z k {v) saturates 
in the IR, in agreement with fig. 3. The proof of mass 
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FIG. 5: The flow of the quantity kd k Log (j!rj$) for 6 = 2 
as a function of k. It saturates in the IR to the value —2, 
showing that ^44 oc A- 



generation in the IR is given by the saturation of the 
curvature V k "{v) provided that Z^(v) stops running. A 

natural quantity to probe this is kdkLog which 

should go to —2 in the case of a finite dimensionful cur- 
vature. As expected, the figure 5 shows that we have 
been able to stabilize a finite mass gap never seen be- 
fore in other renormalization procedures. As we saw, 
its obtention heavily relies not only on truncating of the 
flow equation (2) non-perturbatively in a consistent man- 
ner (3) to access the strong-coupling regime, but also on 
not using the relevance classification that is UV i.e. the 
Fourier expansion of F^ . It illustrates nicely the fact that 
the relevance classification can be totally changed in the 
IR even for a massive theory. We also emphasize that 
the topological excitations responsible for the mass gap 
were never explicitly encoded in this study, showing the 
power of fRG even away from criticality of a Kostcrlitz- 
Thoulcss transition [131. 




flow. We saw on fig. 2 that the flow in the concave re- 
gions lives near the limit (7). We have not been able to 
stabilize the trajectory as far as one wants in the IR as in 
[10] but one can compute the /3-function of a given by (7) 
which corresponds to the evolution of the distance of the 
system to the instability. We found that this /3-function 
is actually suppressed when the mass gap is built, sug- 
gesting the stabilization above the unstable line a = 0, 
see fig. 6. Intuitively, this probable IR stability can be 
understood quite simply. In the IR, the regulator mostly 
behaves as a mass term so that the EEA actually solves 
a serie of massive sGM whose mass parameter is progres- 
sively decreased. As it is well known, the massive sGM 
has the same structure of flows as the sGM for scales 
above its mass so that the interpolation between these 
models is smooth, corresponding to a cross-over located 
at lower and lower energies. 

All the results presented here are stable against a 
change of regulator (changing b or taking an exponential 
suppression) or the use of the propertime scheme [14] 

The author wishes to thank J. -P. Blaizot, B. Delamotte 
and J. Polonyi for many fruitful discussions. 



FIG. 6: The /3-function of a as defined in (7) as a function of 
the scale k. 



We conclude with the problem of the stability of the IR 
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